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ABSTRACT 



Context. Nanoflares are small impulsive bursts of energy that blend with and possibly make up much of the solar background emission. 
Determining their frequency and energy input is central to understanding the heating of the solar corona. One method is to extrapolate 
the energy frequency distribution of larger individually observed flares to lower energies. Only if the power law exponent is greater 
than 2 is it considered possible that nanoflares contribute significantly to the energy input. 

Aims. Time sequences of ultraviolet line radiances observed in the corona of an active region are modelled with the aim of determining 
the power law exponent of the nanoflare energy distribution. 

Methods. A simple nanoflare model based on three key parameters (the flare rate, the flare duration, and the power law exponent of 
the flare energy frequency distribution) is used to simulate emission line radiances from the ions Fexix, Caxm, and Sim, observed by 
SUMER in the corona of an active region as it rotates around the east limb of the Sun. Light curve pattern recognition by an Artificial 
Neural Network (ANN) scheme is used to determine the values. 

Results. The power law exponents, a « 2.8, 2.8, and 2.6 are obtained for Fexix, Caxm, and Sim respectively. 
Conclusions. The light curve simulations imply a power law exponent greater than the critical value of 2 for all ion species. This 
implies that if the energy of flare-like events is extrapolated to low energies, nanoflares could provide a significant contribution to the 
heating of active region coronae. 
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1. introduction 

Heating the corona b y the dissip ation of current sheets 
was first suggested by iGoldl (Il964l) and later devel oped to 
form the ba sis of t he nanoflare heating model by iLevind 
d 1974b and iParked (I1983L Il988h . The idea is that cur- 
rent sheets arise spontaneously in coronal magnetic fields 
that are braided and twisted by random photospheric foot- 
point motions. These current sheets dissipate in many small- 
scale reconnection events, heating and accelerating plasma 
in the coronal loops. In the corona, they would give rise 
to multiple unr esolvable loop strands w i th specific observ- 
able signatures dZirker & Cleveland! 1 19941; I Warren et al.|j 2QQ2: 
ICareill & Klimchukl 12004 iPatsourakos & Klimchukl l2005h . 
Recently lAschwandenl (12008b found evidence against such 
multi-temperature strands in TRACE coronal images. He con- 
cludes that nanoflare heating is only possible if it occurs in the 
chromosphere/transition region where heating across magnetic 
field lines can produce the isothermal loops seen in the corona. 
Irrespective of where the nanoflare energy input sites are, a key 
question is whether the energy of nanoflares is sufficient to heat 
the corona or not. Most of the individual nanoflares would be too 
small to detect and the majority would be small fluctuations on 
the overall background. That background could be produced by 
the blending of many small events. 

The approach taken to estimate their contribution has been to 
extrapolate the energy frequency distribution of detectable flare- 



like events. The energy frequency distribution of larger flares 
tends to follow a power law distribution 



dN 
~dE 



E~ a , 



(1) 



where dN is the number of flares per ene rgy interval dE . The 
energy of small flares dominates if a > 2 (IHudsonlll99l|) . This 
is therefore a critical parameter for the nanoflare heating model. 
The standard method to determine a is to evaluate the energy of 
many flares in a series of observations and then plot their fre- 
quency in bins of energy dE. The majority of analyses based 
on this type of event counting deduce a ~ 1. 7 dLin et al.| [l984: 
IShimizul [19951: lAschwanden & Parnelll 12002b . a value smaller 
than the criti cal 2. These re sults may, however, be misleading. 
For example, [Parne g d2004b demonstrated that one can obtain a 
ranging from 1.5 to 2.6 for the same data set using different but 
still reasonable sets of assumptions for the analyses. 

Here we take an alternative approach and model ultraviolet 
(UV) radiances observed by the Solar Ultraviolet Mea s ureme nts 
of Emitted Radiation (SUMER; IWilhelm etaflll995l 1 1997b in 
an active region corona, assuming that the radiance fluctuations 
and the nearly constant 'bac kground' emission are caused b y 
small-scale stochastic flaring (IPauluhn & Solankil [20041 120071) . 
The model has been app lied successfully to UV ra diance fluctu- 
ations in the quiet Sun (IPauluhn & Sol anki 2007]). The method 
compares light curves generated assuming random flaring with 
a power law frequency distribution to the light curves of an ob- 
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Fig. 1. EIT 195 A images of the observed active region at two 
times, showing the position of the SUMER slit, indicated by the 
vertical line. 
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Fig. 2. Time series of line radiance along the SUMER spectrom- 
eter slit for the period 16 Sep 19:00 UT to 17 Sep 7:36 UT. The 
distance in pixels along the SUMER slit is shown on the vertical 



served emission line. It has the advantage that it takes into ac- 
count without bias weak, blended micro- and nanoflares that pro- 
duce a nearly continuous background. 

Here we apply this technique to off-limb time series 
recorded by SUMER. The three lines modelled, Fe xix X 1 1 1 8.07 
(6.3 MK), Caxm A1133.76 (2.2 MK) and Sim ^ 1113.23 
(0.06 MK), cover two decades of formation temperature from 
the lower transition region to the hotter gas in the corona. 

The analysis described here uses Artificial Neural Networks 
(ANNs) to find the optimum match to the three parameters 
of the model. The main advantage of this method over pre- 
yious analyses based on the radiance di stribution function 
dPauluhn & Solankill2007l: ISafari et al.ll2QQ7h is that we are able 
to obtain quantitative values for all parameters, including a. 
Another advantage of the ANN method is that it concentrates 
on the number and shape of the emission peaks along the light 
curves with little weight on the low r adiance pixels, which was 
a problem with the Saf ari et aD (120071) analysis. 



2. SUMER data and analysis 

The observed active region (AR 1967) is shown in Fig. [II This 
is the region and data set discussed in IWang et al.l (120061) . The 
SUMER 300" x 4" slit was placed, as shown, at a fixed posi- 
tion above the limb. Observations with a cadence of 90 s in six 
spectral lines, Fexix X 1118.07 (6.3 MK), Caxv X 1098.48 and 
,1555.38 (3.5 MK), Caxm X 1133.76 (2.2 MK), Nevi ,1558.62 
(0.3 MK) and Sim ^ 1113.23 (0.06 MK) were transmitted, 
for periods of 12.6 hours followed by a full spectrum scan 
(800-1600 A) of 3.4 hours. A typical time sequence in any one 
line consists of 500 exposures. The three strongest lines, Fexix, 
Caxm, and Sim, are analysed here. Images of their radiance 
along the slit are shown in Fig. [2 for a typical 12.6 hour period. 
Distinct events can be seen in Fe xix, but only the very strongest 
make an impres sion on the bright acti ve region Ca xm emission 
when they cool ( Innes & Wang|l2004b . Si in is seen close to the 
limb and appears to be generated by small surge-like ejections. 
Our results are based on three such time series, taken over the 
days 16-18 September 2000. 

The emission along each row was very noisy at several po- 
sitions. To improve the signal-to-noise but at the same time not 
to lose individual structures, the light curves were obtained by 
first averaging SUMER data over five spatial pixels (5") along 
the slit. Only the light curves with all 500 data points above a 
chosen threshold were selected for analysis. We did not want to 



base the threshold on an absolute intensity because this would 
have biased the input data against low background. So for the 
Fe xix and Ca xm, the threshold was set such that the ratio of the 
Ne vi to Fe xix intensity was less than 0.5. This ensures that only 
light curves from the central part of the active region were taken. 
Most of the Si m emission was concentrated near the limb to the 
south (Fig. [2]). The Sim selection was based on the local scatter 
in the second moment of the line, the line width. If the stan- 
dard deviation of the line width was greater than 1.0 over a local 
5x5 space-time block, then the central data point and associated 
light curve were excluded from the analysis. This resulted in 35 
test light curves for both Fe xix and Ca xm and 1 1 test curves 
for Si m. Before being fed to the neural network, all light curves 
were normalized to their maximum. 

3. Model 

3.1. Applicability 

The emission in the active region corona is assumed to be caused 
by many random flares with flare radiances following a power 
law frequency distribution. Flares with a power law frequency 
distribution, a, in radiance are assumed to erupt with a fre- 
quency, pf, and have a flare duration r = r r + r^, where r r is 
the rise time and the decay time. We assume r r /rj = 0.5. The 
other free parameter in the model is the ratio of the maximum to 
minimum flare energy which is set to E m3iX / E m { n = 50. 

For a large number of independent random flares, the distri- 
bution of normalized radiances, / = I /I where I is the radiance , 
is lognormal with shape parameter cr dPauluhn & Solankill2007l) : 



f(J) 



1 



crJ^Zn 



exp 



dog J) 2 

2(T 2 



(2) 



cr is inversely proportional to \Jrpf (Pauluhn & Solanki 2007), 
with a slight a dependence ( Safari et al .1120071) . A small shape 
parameter (cr < 0.3) indicates a symmetric distribution due to 
high background emission caused by either a long duration time, 
r, or a high flare frequency, p/. The radiance distributions of the 
three lines of Fe xix, Ca xm and Si m and their lognormal fits 
are shown in Fig. [3] This gives us confidence that the stochastic 
flare model is applicable. It is interesting to note that both the 
Fe xix and the Si m lines have the same shape parameter. 

Light curves for the stochastic flare model are shown for a = 
1.6 and a = 2.4, and two combinations of rpf = 2 in Fig. [4] The 
light curves are visibly different, although they all have shape 
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Fig. 3. Distribution functions of SUMER data in the active re- 
gion corona (solid lines) and best fit lognormal functions (dashed 
lines). The radiances are normalized to their median and their 
distributions to the number of data points. 
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Fig. 4. Light curves for flare models run with different r, pf and 
a parameters. All light curves have rpf = 2. 



parameter cr « 0.6. The effect of a on the light curve is seen in 
the ratio of strong to weak flares. The left-hand light curves have 
more large flares because they have a smaller exponent. Picking 
up these pattern changes is the strength of the ANN method. 



3.2. Neural networks and parameter estimation 

We applied the ANN method to probe the unknown parameters 
(power law exponent, a, duration time, r, and flare rate, pf) of 
the three lines. ANNs have become a popular tool in almost ev- 
ery field of science. In recent years, ANNs have been widely 
used in astr onomy for applications such as star/galax y dis- 
crimination, (lAndreon et al.ll2000t ICortiglioni et al.ll200ll). mor- 
pholo g ical classificatio n of galaxies. (iStorrie-Lombardi et al.l 
1992; Balletal. 2004) and spectral classificati o n of stars 
(Ivon Hippel et al.lll994i iBazarghan & Guptal l2QQ8t iBazarghanl 
120081) . 



W e emp loy probabilistic neural networks (PNNs ISpechtl 
119881 ll99Qh . The PNN learns to approximate the probability 
density function of the training samples. It uses a supervised 
training set to develop distribution functions within a pattern 
layer. These functions in the recall mode are used to estimate 
the likelihood of an input feature vector being part of a learned 
category or class. 

An example of a PNN is shown in Fig. [5] This network has 
four layers. The network contains an input layer which has as 
many elements as there are separable parameters needed to de- 
scribe the objects to be classified. It has a pattern layer, which 
organizes the training set such that each input vector is repre- 
sented by an individual processing element. The third layer is 



Fig. 5. Schematic of a typical probabilistic neural network. 



the summation layer, which has as many processing elements 
as there are classes to be recognized. Each element in this layer 
combines via processing elements with the pattern layer which 
relates to the same class and prepares that category for output. 
Finally, there is the output layer that corresponds to the summa- 
tion unit with the maximum output. 

For the identification of SUMER light curves, the input vec- 
tor, X = (x\,X2, ...,x n ), is the light curve with 500 data points 
(n = 500). The network is first trained to classify light curves 
corresponding to all the possible combinations of a, r, and pj. 
For this we synthe tically generate light curve s with the nanoflare 
code described in iPauluhn & Solankil (120071) . We generate one 
light curve for each combination of the parameters: 

- the power law exponent spanning 1.5 < a < 3.2 in steps of 
0.1, 

- the duration time spanning 1.5 < r < 52 in steps of 1, 

- the flare rate spanning 0.1 < p/ < 0.9 in steps of 0.1 with 
additional values at 0.05 and 0.95. 

This gives a set of 6930 pattern groups (k = 6930), one group 
for each combination of a, r, and pf. E ach pattern group, k , is 
characterized by Nk Gaussian functions (ISpechtll 1988lll990h . 

When a SUMER light curve of an unknown classification is 
fed to the network, the summation layer of the network computes 
the probability functions Sk of each class. Finally at the output 
layer we have C, the value with the highest probability. 



4. Results and conclusions 

4.1. Results 

In the present work, PNN is used as a tool to extract the three 
flare model parameters required to reproduce the SUMER light 
curves. All 35 Fexix and Caxm, and 11 Sim SUMER light 
curves from the three days of observations were fed individu- 
ally into the neural network and the parameters were obtained 
for each light curve separately. The final PNN outputs are shown 
in Table [T] The bold numbers are the statistically maximum oc- 
currence for each parameter. For example for Fexix, a = 2.8 is 
found in more than 70% of the light curves. The minimum and 
maximum values, given on the left and right, indicate the scatter 
in the light curve parameters. 

In each line there is 20% scatter in a, and 50% scatter in 
r. The range of pf values for Fe xix and Si m is much broader, 
suggesting that events producing emission in these temperature 
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Table 1. The SUMER spectral lines and the parameter values 
given by PNN. 



SUMER spectra 




PNN outputs 1 




lines 


a 


T 


Pf 


Fe xix 


2.5 2.8 3.1 


9 14 20 


0.1 0.2 0.7 


Ca xm 


2.6 2.8 3.0 


41 45 45 


0.8 0.9 0.9 


Sim 


2.4 2.6 2.8 


5 9 12 


0.2 0.3 0.9 



1 The most frequent values are given in bold, and the min- 
imum and maximum on the left and right, r is per expo- 
sure time (90 s) and p f is per exposure time per 5"x4" 
spatial element. 



ranges do not have the same rate everywhere but are seen in ir- 
regular bursts. We also note that the value of rpf is roughly the 
same for both Fe xix and Si m, as suggested by their shape pa- 
rameter (Fig. [3]). The Caxm light curves are all matched with 
a high value of pf, consistent with the idea that the 1 MK ac- 
tive region corona requires almost continuous flaring. The four 
times higher rate for Ca xm than Fe xix suggests that most of the 
Ca xm emission is produced by heating events below the Fe xix 
formation temperature (6.6 MK). 

Example light curves obtained using these parameters are 
compared with the observed ones in Fig. [6) Both the Sim and 
Caxm simulations look remarkably similar to their observed 
light curves. The background radiance of the Fexix light curve 
is about a factor of 2 too low. The Fexix light curves had a pf 
ranging from 0.7 to 0.1, so we suspect that in this case the pf 
value is slightly too low. Also for Fe xix, the ratio r r /rj deduced 
from the data is smaller than the fixed value 0.5 used here. This 
may influence the accuracy of the method. 

The sensitivity of the PNN output depends on the training 
set. During the training session, the network must see all possible 
patterns that it is supposed to classify in the testing session. With 
500 simulated light curves in the training set, PNN was not able 
to converge for several of the SUMER light curves. When we 
increased the number of simulated light curves to 6930, we were 
able to obtain unique parameters for all observed light curves. 

4.2. Conclusions 

The concept that the solar corona may be heated by numerous, 
randomly distributed, small flare-like events called nanoflares is 
considered by comparing simulated and observed emission line 
light curves. The difference between this and previous methods 
is the fully automated modelling of the light curve structure. 
There is no human decision required for background/event cut- 
off levels or best fit parameters. 

The result is power law flare energy frequency exponents 
greater than 2.5 for all three emission lines considered, Sim, 
Caxm and Fexix. This is consistent with the corona being 
heated mainly by nanoflares, and demonstrates the importance 
of nanoflare 'background' emission in determining the power 
law exponents. The parameter with highest uncertainty or largest 
scatter is the flare rate, especially for the lines formed at transi- 
tion region and hot flare temperatures. Coronal plasma at these 
temperatures is produced sporadically and is associated with 
more specific coronal and chromospheric loop structures than 
the general active region corona, so the scatter is to be expected. 

The next step will be to determine the actual flare energies 
producing the nanoflare emission. This is a much more compli- 
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Fig. 6. Samples of the radiance time series: left panel: SUMER 
data, and right panel: simulation data obtained with the parame- 
ters given in Table 1 . 

cated exercise because the modelled light curves are observed in 
the corona w hich may be heated by events occurring lower in the 
atmosphere (lAschwa nden 20Q8|), so that it requires a model for 
the energy transfer to the observation position. 
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